arXiv:gr-qc/9808011vl 4 Aug 1998 


Phase diagram of the mean field model of 
simplicial gravity 

P. Bialas a * Z. Burda 6 and D. Johnston 0 


“Fakultat fur Physik, Universitat Bielefeld, 33501 Bielefeld Germany 
^Institute of Physics, Jagellonian University, 30-059 Krakow, Poland 
c Dept, of Mathematics, Heriot-Watt University, EH14 4AS Edinburgh, Scotland 


Abstract 

We discuss the phase diagram of the balls in boxes model, with 
a varying number of boxes. The model can be regarded as a mean- 
field model of simplicial gravity. We analyse in detail the case of 
weights of the form p(q) = q ~^, which correspond to the measure term 
introduced in the simplicial quantum gravity simulations. The system 
has two phases : elongated (fluid) and crumpled. For j3 £ (2, oo) the 
transition between these two phases is first order, while for /? £ (1, 2] it 
is continuous. The transition becomes softer when f3 approaches unity 
and eventually disappears at /3 = 1. We then generalise the discussion 
to an arbitrary set of weights. Finally, we show that if one introduces 
an additional kinematic bound on the average density of balls per box 
then a new condensed phase appears in the phase diagram. It bears 
some similarity to the crinkled phase of simplicial gravity discussed 
recently in models of gravity interacting with matter fields. 


The mean field description of simplicial gravity introduced in Jl], || |J was 
analysed in the series of papers 1 0 0. As discussed there, the mean 
field approximation explains many facts observed in numerical experiments 
of simplicial gravity, such as the appearance of singular vertices and the 
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mother universe, and the discontinuity of the phase transition. Recent work 
00 shows that adding vector matter fields to simplicial gravity effectively 
amounts to renormalising the a parameter in the measure term q a where q is 
the vertex order. The resulting phase diagram resembles the phase diagram 
of the mean-field model which we discuss in detail in this letter. 

The model is given by the partition function || : 

Mmax 

Z(N,k, p) = Y e * M E J>(91) 

M= 1 qi,—,qM 

Mmax 

= Y e* M z(N,M,l3) 

M= 1 

It describes an ensemble of N balls distributed in a varying number of boxes, 
M. The system has at least one box and at most M max boxes. The partitions 
of balls are weighted by the product of one-box weights p(q), where q is 
the number of balls in the box. Here we consider one parameter family of 
weights : 

p(q) = q~ p for g = i,2,.... (3) 

Note that for these weights the minimal number of balls in a box is q m i n = 1. 
If there are no further constraints this implies that the maximal number of 
boxes is equal to the number of balls M max = N. 

For large N the partition functions z (...) and Z(. ..) are expected to 
behave as : 

z(N, M, (3) ~ N v e N ^ r, P > with r = — , (4) 

and 


• • -p(qM)Sqi 4 - \-q M , N ( 1 ) 

( 2 ) 


Z(N, k,/3) ~ N Y e N ^ K ^ 


(5) 


where f(r,/3 ) and 0(/c, /3) are appropriate thermodynamic potential densities 
depending on the intensive quantities characteristic for the (r, (3) and (k, /3) 
ensembles. The power law corrections for the partition functions with the 
exponents y, Y disappear in the thermodynamic limit. In our case they 
are the leading corrections. In general, however, there might be stronger 
corrections, for example of the type expcIV“, with a < 1. 

Substituting the summation in (|2|) by an integration over a continuous 
variable, we can relate the partition functions by the transform : 


N Y e N( P( K >P) — _/\r y+1 



( 6 ) 
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The additional factor N on the right hand side of (p[) comes from the inte¬ 
gration measure Ndr. In the thermodynamic limit, this gives : 

= f{r*,P) + r*Ac. (7) 

where r* is maximum of the integrand in (|G|). The value of r* is given by the 
equation : 


- k = d r f(r„(3) (8) 

if the maximum of the integrand lies inside the interval (0,1) and is equal 
r* = 0 or r* = 1 otherwise. Additionally we see that in the case when 
the maximum lies inside the interval (0,1) the integration over r around 
the saddle point introduces the additional factor N~ 1 ^ 2 . This leads to the 
following relation between the exponents y and Y : 

Y = V+\- (9) 

Otherwise, the maximum lies at either end of the interval and the integration 

over r introduces a factor iV _1 , thus 

Y = y. (10) 

We shall return to the exponents Y, y later, concentrating now on the ther¬ 
modynamic limit. 

In the (k, /3) ensemble one defines the average (r) : 

(M) 1 d K Z(N , k, (5) 

{r) N N Z(N, k, f3) ' 1 ’ 

As we will see this quantity plays the role of the order parameter. In the 
thermodynamic limit the equation ( |TT| ) reads : 

(r) = d K <f>(K, (3) = r*. (12) 

The second equality in the last formula can be treated as an inverse transform 
to (1). It results from the fact that the integrand of (H), which corresponds 
to the distribution of r, is sharply peaked around r* when —» oo and the 
position of the peak also becomes the average of the distribution. 

The task of determining r* relies on finding the maximum of the function 
f(r,/3) + nr. The function f(r,/3) can be found by the saddle point method 
0 : 


/( r,0) = (i,(r) +rK(ii,(r),0) 


(13) 
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where /x*(r) is a value which maximises the right hand side of the above 
expression for the given r. The function K(n,/3) is defined as : 

OO 

= log J2p(q) e ~ fig ■ ( 14 ) 

9=1 

The properties of the function K (/x, /3) are central to the further consider¬ 
ations as they fully determine the behaviour of the model. For the particular 
choice of weights (|3|) : 

OO 

K(/d,p) = log = log g p (fj) (15) 

9=1 


where gp{/i) is a function familiar from Bose-Einstein condensation, which 
has an integral representation : 


1 r°° t^dt 

w) - J o ell+t _ x • 

(16) 

The function g^(fji) is defined on the interval [0, oo), 
of the saddle point equation : 

so /x* is either a solution 

i = ~d„K(v„l3) = 
r gpKn*) 

(17) 

for r > r cr or else //* = 0. The critical value of r is 


( 0 (3< 2 

r cr — \ C(/3) a ^ n , 

l CC9-1) P 

(18) 

where ((/3) is the Ricmann Zeta function which arises when one inserts /x* = 0 
into K(/i, (3) in ([15]). For r < r cr the function g, + rK(n, (3) has no maximum 
in the range (0, oo) and /i* = 0. 

We thus see that the function f(r,/3) as a function of r changes regime 
at r = r cr . For r < r cr ie fi* = 0 the equation (|13|) reduces to the linear 
dependence : 

f(r,P) = -Ka.(P)r 

(19) 

where 


- K cr = logC(^) • 

(20) 


We have added the minus sign in the definition of K cr for convenience to 
adjust to the sign in the saddle point formula (||). The function f(r,/3 ) as a 
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Figure 1: The f(r,/3 ) for various values of (3. The curves for (3 > 2 possess 
a linear part at r — 0. The main difference between the curves for (3 = 0 
and (3 = 1.5 is the former has an infinite derivative d r f at r = 0, whereas 
the latter has a finite derivative. For all the curves the derivative at r = 1 is 
infinite. See figure [^. 

function of r has different qualitative behaviour depending on (3. For (3 > 2 
the function has a linear piece for r < r cr . This piece becomes shorter when (3 
approaches 2 and finally disappears. For 1 < (3 < 2 the function has no linear 
piece, but it does have a finite derivative at r = 0. This distinguishes this 
range from (3 < 1 where the derivative d r f at r = 0 is infinite. Representative 
shapes of the function / for (3 from the three ranges are shown in figure [l]. 
This figure is supplemented by figure |] where the derivative is plotted to 
expose its behaviour at r = 0. With this picture in mind one can see the 
main feature of the solution r* of the saddle point equation (|]). Namely, for 
(3 < 1, for each k it has always a solution in the range (0,1). For 1 < (3 < 2, 
r* lies in the range (0,1) so long as n > K cr , whereas r* = 0 for n < K cr . 
Finally, for (3 > 2, the solution lies in the range (r cr , 1) for k > n cr and then 
jumps to r* = 0 for n < n cr . This is illustrated in figure [|. The discontinuity 
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Figure 2: The derivative d r f(r,f3 ) of the curves from the figure [T]. All go 
to infinity at r = 1. Depending on the value of /3 they drastically change 
behaviour at r — 0 : for (3 = 0 they go to infinity, for all other are finite. For 
/3 = 2.25 and [3 = 3.0 they are constant in some range of r. 


at k — K cr (pOD has the height r cr given by (|l^). The average (r) = r* is 
an order parameter for the model. A simple consequence of equation ( |TE|) 
is that the transition becomes continuous when (3 — 2 since r cr = 0. In 
fact, the transition stays continuous as long as (3 is in the range (1,2]. In 
this range the r cr = 0 and n cr is finite and given by (|20|). As (3 approaches 
unity, K cr —> — oo and the transition disappears. This process is illustrated 
in figure [|. For (3 < 1 the values of the derivative d r f(r,/3) span the interval 
(—oo, oo) and the saddle point equation (||) always has a nonzero solution 
r* for the whole range of k. The model therefore has only one phase and no 
phase transition. To illustrate this, consider the solution when (3 = 0 : 


r* = 


1 + e" 


( 21 ) 


which is a smooth function spanning the range (0,1) without any singularity. 
The phase diagram in the (k, (3) plane is drawn in figure H. The solid 
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Figure 3: The value of v*(k) for various values of (3. The portions of plots 
lying on the x axis (r* = 0) are not plotted. 


line is a discontinuous phase transition with a nonzero gap r cr (|T8|) while the 
dashed line is a continuous transition. 

Coming back to the range (1,2) of /3 one can determine the singularity 
type of the thermodynamical functional <J)(k, /3) at the transition ie when 
k —> n cr . For small fi : 

o° _ k 

gM = ni- 0 )/-' + Y,<; IP-Q-n-- (22) 

k =0 

Inserting this into the equation ([TtI) one obtains for small /i : 

r ~ fi 2 ' 13 + ... , (23) 

or, after inverting it for fi : 

^ ~ . (24) 


In this range of (3 r cr = 0 and so /i* is given by the equation : 

K = -log^(/x*) (25) 
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Figure 4: The phase diagram. The continuous line denotes a first-order phase 
transition. The dashed line denotes a continous transition. The dotted line 
shows the transition line that appears after addition of the artificial cut-off. 


For small /i* we can write : 

K- K cr = - log gpifj,*) + logC(/3) ~ // _1 ~ r (d-h/(2-/3) _ (26) 


Inverting this with respect to r we eventually get : 


(r) = r* ~ (k - K cr ) x , 


where 


x = 


2-/9 

/9-1 


(27) 


The exponent x varies between oo and 0 as f3 grows from one to two. In the 
limiting case (3 — > 1, the exponent x — > cxo which means that the singularity 
becomes arbitrarily soft, which perfectly matches the fact that for (3 < 1 there 
is no singularity. At the other end of the interval, when (3 —> 2, the exponent 
x —» 0, which means that (r) becomes an arbitrarily steep function at n cr 
which again agrees with the fact that for [3 > 2 the function is discontinuous 
at n cr . 

What we have discussed so far is an extension and refinement of the 
analysis of the model presented in [J3J. A few comments are in order. The 






analysis can be generalised for an arbitrary set of weights. In general, the 
transition takes place when the argument e _M of the series fll4l) approaches the 
radius of convergence of the series. Denote the radius by e A . In our example 
the radius was 1 ie A = 0, but in general A may be any finite number. The 
condition for the existence of the transition is that the generating function 
K (fT4l) is finite at the radius of convergence of the series : 

K{ A, f3)< oo . (28) 

The family of weights may be parametrised by many parameters. In the last 
formula, the parameters are represented by a single symbol f3. The equation : 

-K cr = K(A,/3), (29) 

is the equation of the critical line between the phases of the (k, /3) ensemble. 
Moreover, as one can see by repeating the arguments previously used for the 
particular weights ([|), the condition for the transition to be discontinuous 
is : 


dpK(A, (3) < oo 


(30) 


since the inverse of the derivative d lt K corresponds to the latent heat. The 
transition is continuous for /3’s for which the latent heat vanishes and condi¬ 
tion (^8|) is fulfilled. 

The position of the critical line on the phase diagram can be easily 
changed by a slight modification of weights. For example, a simple mul¬ 
tiplication by a constant C : p(q) —> Cp(q) corresponds to a horizontal shift 
of the critical line in the phase diagram (Fig. }4]). In general a slight modifi¬ 
cation of weights can also change the (3 position of critical line and thus the 
critical exponent (|27|). In this respect this model is “non-universal” (see also 

!)• 

So far we have discussed the behaviour of the model in the thermodynamic 
limit. We now determine the exponents y (||) and Y (^). It is convenient to 
define the grand canonical partition function : 

OO 

Z(/i, «,/?)=£ e~ N ^Z(N , k, (3) ~ {p - ^o)- 1 '"', (31) 

N =1 


whose singularity type depends on the exponent Y. The critical value of the 
chemical potential is /i 0 = 0 (k, (3) is zero in the crumpled and nonzero in the 
fluid phase. The value of the exponent Y depends also on the position (/c, (3) 
in the phase diagram. From (Q) one obtains the following expression : 


Z{l*,K,(3) 


1 - e K gp(p) 


(32) 
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from which one can extract the singularity. In the fluid phase, the singularity 
comes from the residual denominator and it is : 

~ (/i-/i 0 ) _1 , (33) 

from which one concludes that Y = 0. 

For (k, (3) in the crumpled phase, the singularity comes from the numer¬ 
ator of the fraction (^) which itself has a power like singularity when /i 
approaches /i 0 = 0. Thus 


~ fJ,' 


P -1 


(34) 


and hence Y = — (3. In other words, the partition function inherits the 
singularity from the weights or more precisely from the generating function 
K(n,(3). This effect comes from the dominance of the singular box contribu¬ 
tion over the negligible entropy. In general the singularity need not be power 
like. For example the weights : p(q) = exp (—aq + bq 1 ^ 2 ), where a and b are 
positive constants, would lead to an essential singularity. In this case N 5 
would not be the dominant correction to the partition function. 

Let us now determine the singularity type on the critical line. The sin¬ 


gularity comes from the denominator in the expression (|32|) . It changes at 
/3 = 2. For (3 > 2 the expansion of the denominator in small /i begins with 
the linear term. Because the numerator approaches a non-vanishing constant 
for /i —> 0, we have : 


Z{^k,(3) ~ [i 1 , (35) 

and hence Y — 0 as in the fluid phase. For 1 < (3 < 2 the expansion of the 
denominator begins with the singular term p 131 (|T6|) , thus 

Z{h,k,P) ~ fi 1 - 0 , (36) 


and Y = —2 + f3. 

So far we have tacitly set q m in = 1, which meant that each box contained 
at least one ball. As a result, the system with N balls could have M boxes 
at most. Hence the maximal value of r is 1. If one instead chooses qmin to 
be larger than one then the maximal r is : 

Tynax 1 / qmin 

In this case, the analysis of the model proceeds exactly as before, except that 
in place of 1 for the upper limit one sets r max = \/q m in■ In particular, the 
range of integration in equation (bj) is now (0 ,r max ) and the derivative d r f 
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is infinite at r max . The average (r) approaches asymptotically r max when 
k —> —oo. So the only essential difference to the previous analysis is the 
range of r. In particular, the phase structure of the model is unchanged. We 
shall call r max the natural limit for r. 

The phase structure of the model can be made more complex by intro¬ 
ducing a slight modification. Let us impose a kinematic upper limit on r. 
We shall call this the artificial limit and denote it by r art . The artificial limit 
obviously has to lie below the natural one, r art < r max . Now the partition 
function becomes : 

T art 

J^ Y e N <j>(K,,l3) _ jyi/+l J ^ re N(f(r,/3)+rK) (gg) 

0 


The main difference between the natural and the artificial limit is that the 
derivative d r f is infinite at the former and finite at the latter. In the case of 
a model with the natural limit, this limit is never reached by the average (r), 
or at worst reached only asymptotically for k —> — oo, whereas the artificial 
limit can be reached for finite k when the term tk dominates over f(r,/3 ) 


in the integrand d38|). If we fix (3 and change k from large positive values 
towards large negative, the system is first in the condensed phase (r) = 0, 
then at the critical value (r) jumps to r cr and later smoothly rises until it 
reaches r art where it sticks (see figure [5]) . Such a situation is very similar to 
the previous standard one. 

There exists, however, the possibility of quite different behaviour. For 
/?>/?* where /?* is defined by the equation : 


^art 


C(A) 

C(ft -1) 


(39) 


the value of r cr is larger than r art . In this case, at the transition from the 
condensed phase (r) jumps from zero directly to r art . It cannot jump to 
r cr because this lies above r art and is beyond the allowed kinematic range 
(figure 0). This gives a new phase which is neither crumpled nor fluid. It 
has (r) > 0 but on the other hand there is a singular box in the system. We 
call this phase condensed , as it is analogous to the condensed phase of the 
canonical Z(M,N,/3) ensemble. We have marked the critical line separating 
the condensed phase from fluid one by the dotted line on the phase diagram 
in the figure [|. 

Define the quantity : 


(s) = 



N 


(40) 
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Figure 5: The dependence r* on k for various values of f3 in presence of the 
artificial cut-off 


which is the fraction of all balls which are in the singular box. The quantities 
(r) and (s) play the role of the order parameters. We have : 

crumpled : (r) = 0, (s) = 1, 

fluid : (r) > 0, (s) = 0, (41) 

condensed : (r) > 0 , (s) > 0 . 

The phase diagram is very similar to the one of simplicial gravity with 
the matter fields or the measure term j7]. Q. If one naively applies the mean 
field model to 4d simplicial gravity one would get the following constraint 
for the orders of vertices : 

Qi + Q 2 + • • • + Qn 0 — 5 IV 4 (42) 

where Nq is the total number of vertices, g* the order of vertex i , and IV4 
is the number of 4simplices. We thus have the correspondence IVo M, 

Ko —k and IV4 N. The minimal order of a vertex is q m i n = 5 , so the 
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natural limit r max = 1/5. This means that : 


m 

N. 


= 1 


( 43 ) 


is the natural upper limit for (No) / N\ from the point of view of the mean field 
model. As we know this limit is never reached since the simplicial structure 
imposes the limit (No) = 1/4 which in the language of the mean field model 
corresponds to the artificial upper limit [fTTH . The full 4d theory thus provides 
the structure which we have assumed in the discussion of the mean-field 
model. Of course, the cut off at 1/4 is not as sharp as we have assumed. One 
expects a smooth cut off, which means that the entropy of configurations for 
(Nq) / IV4 approaching 1/4 gradually decreases in comparison with the mean 
field model, so that the upper limit (A'q) /N 4 = 1/4 is not reached at finite 
Ko but rather asymptotically approached for large /to- 

The source of the limit (A r o)/A r 4 = 1/4 in simplicial gravity is purely 
geometrical and results from the fact that qt is the number of 4-simplices. 
A consequence of this is that if we add a point with qi = 5, for example 
by a barycentric subdivision, we increase the numbers q 3 of neighbouring 
vertices|ll|]. In other words, q t 's are not independent and in particular one 
cannot make all of them equal to 5 as the natural limit would require. 

The exponent Y requires some comment. In the theory of random geome¬ 
tries one defines the string susceptibility exponent 7 , which is the counterpart 
of Y. The exponents Y is shifted by a constant with respect to the standard 
definition of 7 : Y — 7 — 3. If one followed these definitions in the fluid 
phase one would obtain 7 = 3 in the (/t,/?) ensemble since, as we have seen, 
Y — 0. However, the mean field model neglects all the subtle effects and 
combinatorial factors coming from topological constraints imposed by the 
local geometry so one should treat this value with some caution. 

Similarly, the crumpled phase of the ball in boxes model has Y = —/?, 
but the following argument suggests that the corresponding value of 7 is not 
that which is naively expected. The balls in boxes model in the particular 
case of the r = 1/2 ensemble can be mapped onto the branched polymer 
model with the topology of sphere p], | 9 j • One can recognise the fluid phase 
of the balls in boxes model as the generic phase of the branched polymer 
model. For the former y = —1/2, whereas for the latter 7 = 1/2, from which 
one can conjecture y = 7 — l[|. Switching back to the varying number of 


boxes ensemble we have Y = y in the crumpled phase from equation (p~0D , 
and 7 = 1 — /?. 


^We have y rather than Y as we are dealing with a fixed number of boxes ensemble in 
the branched polymer analogy. 
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To calculate 7 and hence Y, y in the condensed phase we will use a heuris¬ 
tic argument. We know that Y = y since the phase is again on the kinematic 
border m as for the crumpled phase, but y is now the exponent of the 
fixed r = r art ensemble. As before, in the particular case r = 1/2 we can 
take the value of 7 = 2 — (3 from the branched polymer model (i). if we 
assume that this relation holds independently of r, and that y = 7 — 1 as 
conjectured above for the crumpled phase, we obtain Y — y — 1 — (3 which 
differs by one from the crumpled phase value. 

We think that the main feature of the mean-field solution is that the 
exponent 7 properly reflects an universal behaviour, ie it is independent on 
the location in the phase diagram as long as the system is in the fluid phase. 
This is more important in the context of simplicial gravity than the particular 
predictions for the exponent 7 calculated for our choice of power law weights. 
Mean field theories will, in any case, generically predict critical exponents 
incorrectly. 

In summary, we have seen that the mean field model of simplicial quantum 
gravity can account for many of the features seen in simulations of the full 
model. With the correct choice of ensemble (a varying number of boxes) we 
reproduce the first order nature of the transition for for (3 G (2, 00 ), while for 
f3 G (1,2] it is continuous, disappearing altogether at f3 = 1. We have also 
remarked that taking further account of the local geometrical constraints in 
the model can be implemented as a lowered limit on the region of integration 
in the saddle point equation, which we have denoted as the artificial upper 
limit. With the artificial limit in place we found the modified phase diagram 
of figure [|. which possesses a new condensed phase. This is similar to the 
crinkled phase of full simplicial gravity with matter fields or the (apparently 
equivalent) measure term J7|, ||. 
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